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Abstract 

The CLEO-c experiment has been collecting data at the charm-threshold re- 
gion. A selection of recent results on charmed meson and charmonia decays 
are presented. 



1 Introduction 



The CLEO-c experiment has been taking data at the CESR symmetric e + e~ 
collider at the charm threshold region since 2003. The main goal of the ex- 
periment is to perform high precision measurements of hadronic branching 
fractions, leptonic decay constants, and semileptonic form factors of charmed 
mesons, together with an extensive study of QCD spectroscopy in the charmo- 
nium sector in order to provide rigorous constraints on the strong interaction 
theory, especially Lattice QCD calculations. If the theoretical calculations sur- 
vive these tests they can be used to provide much needed theoretical input to 
extract quark mixing (CKM) matrix elements (e.g. V u b, Vtd and Vt s ), which 
remain limited by complications caused by strong interaction dynamics. 

The selected topics discussed here are the measurement of the absolute 
branching fraction of Cabibbo- favored hadronic D°, D + , and D s decays; mea- 
surement of the leptonic decays, — > l + v, and decay constants /d (s) ; mea- 
surement of the Do mass; and a study of three-body hadronic decays of Xcj- 

Charged and neutral D mesons are produced at the ip(3770) which pre- 
dominantly decays to D + D~ and D°D° with a total cross section of about 
7 nb. D s mesons are created at around E cm = 4170 MeV, where their pro- 
duction is dominated by e + e _ D*^Df with a cross section about 0.9 nb 

1). The l 3 Pj (J = 0, 1, 2) charmonium states are produced in radiative "0(25) 
decays with a branching fraction of 9% to each. The main advantage of the 
CLEO experiment compared to B factories and fixed target experiments is the 
very clean experimental environment with low multiplicity final states, which 
arises from running at or slightly above production thresholds. Background is 
further reduced in e + e~ — > ^(3770) — > DD and e + e~ — > D*D S data by fully 
reconstructing (tagging) one of the D( s ) decaying into a hadronic final state. 

2 Absolute D°,D + ,D S hadronic branching fractions 

Precise knowledge of the absolute hadronic branching fractions of the D°, D + , 
D s mesons is important because they are used to normalize the decays of other 
charmed mesons and mesons. 

CLEO measures the absolute branching fraction of three D°, six D + , 
and six £>+ Cabibbo-favored hadronic decays using single tag and double tag 

events following a technique pioneered by the MARK-III Collaboration © . In 
single tag events only one of the -D( s ) or is reconstructed in a specific final 
state, while in double tag events both and mesons are reconstructed 
in one of the hadronic final states. The single and double tag yield can be 
expressed as rii = NuuBiCi and riij — NjjjjBiBjeij, where Ndd is the number 
of D°D°, D + D~, or D+D~ events produced; Bi is the branching fraction of 
decay mode i; ei and ey are the single and double tag efficiencies. Then the 



absolute branching fractions can be obtained from the double and single tag 
ratios and efficiencies as 

Li = (1) 

rij 

Since ejj rs the branching fraction is nearly independent of the efficiency 
of the tagging mode, and many systematic uncertainties cancel in the ratio. 

There is a difference in the kinematics of D and D s mesons. The D and 
D mesons produced in e + e~ — > -0(3770) — > process have the same well 
defined energy (and momentum) in the center of mass frame of the colliding 
e + e~ beams (Ed = -Ebcam)- In contrast, a pair of D s mesons is produced in 
e+e" -> D^DJ followed by the decay D*± -> jDf (96%) or n°Df (4%). 
Therefore, the D s produced directly has a well defined energy and momentum 
in the e + e~ center of mass frame, while the secondary D s from the D* decay 
has a much broader momentum distribution around the same value. This 
difference in kinematics leads to a slightly different selection strategy of DD 
and DfDJ events. 
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Figure 1: Beam- constrained mass distribution of D (D) candidates in double 
tag events summed over all decay modes. 

In order to identify (tag) the D mesons, we use AE = Er> — -Ebeam and 
the beam-constrained mass, M^ c = y^beam — (pd) 2 , where Eu and po are the 
energy and three-momentum of the reconstructed D meson candidate. Substi- 
tuting the beam energy for Ed improves the mass resolution of D candidates 
by an order of magnitude, to about 2 MeV/c 2 . AE peaks around zero and M\, c 
peaks at the nominal D mass. We require AE to be consistent with zero within 
3 standard deviations, and extract the number of single and double tags from a 
fit to the one-dimensional and two-dimensional Afbc distributions, respectively. 
Fig. Q] illustrates the beam-constrained mass distribution for double tag events 



summed over all decay modes. In 281 pb 1 data, we reconstruct 230, 000 sin- 
gle tag and 13, 575 ± 120 double tag D°D° events, and 167, 000 single tag and 
8, 867 ± 97 double tag D + D~ events. 

The D° and D + branching fractions are determined from a simultaneous 
least squares (x 2 ) fit to all D° and D + single and double tag yields. The fit 
properly takes into account correlations among all statistical and systematic 
uncertainties. The preliminary branching fractions based on 281 pb -1 data are 
listed in Tableland compared to the 2004 PDG averages ®, which does not 
include our earlier results based on 56 pb -1 data ®, in Fig. [2] 

Table 1: Preliminary D° and D + branching fractions with statistical and sys- 
tematic uncertainties. 



Decay 


B(%) 




K~TT+ 




3.88 ±0.04 ±0.09 


D° -> 


K-JT+ 


7T° 


14.6 ±0.1 ±0.4 


D° -> 


K-JT+ 


7T~7r + 


8.3 ±0.1 ±0.3 






-7T+ 


9.2 ±0.1 ±0.3 


D+ -» 




-TT+1T 


6.0 ±0.1 ±0.2 


D+ -> 


K° s ir+ 




1.55 ±0.02 ±0.05 


D+ -> 


K° s ir+ 


7T° 


7.2 ±0.1 ±0.3 


D+ -> 


K° s tt+ 




3.13 ±0.05 ±0.14 


D+ -> 


K-K 


+ 7T+ 


0.93 ±0.02 ±0.03 



In order to select D^DJ events, we use the beam-constrained mass 
(M bc = v^boam _ (pd s ) 2 ) and the invariant mass (M(D S ) = \Je 2 Ds - (pd s ) 2 ) 

of the D 8 (or D s ) candidate and ignore the 7 or ir° resulting from the D* decay. 
The beam-constrained mass is used as a proxy for the momentum of the D s 
candidates (see Fig. [3]). We apply a cut on Mbc that selects all of the directly- 
produced D s and, depending on the decay mode, all or half of the secondary 
D s . Then the invariant mass of the D s candidate is used as a primary analysis 
variable to extract the number of tags. Single tag yields are obtained from 
fitting the one dimensional M(D S ) distributions, while double tag yields are 
determined by counting events in the signal regions in the M(D+) vs. M(D~) 
plane and subtracting backgrounds estimated from sideband regions. 

For this analysis, we use a binned likelihood hybrid fitter which utilizes 
Gaussian statistics for single tag modes and Poisson statistics for double tag 
modes, since the least squares x 2 fitter used for the D branching fraction mea- 
surement is not appropriate for the small signals and backgrounds in the D s 
double tag samples. The preliminary branching fractions based on 195 pb -1 
data are summarized in Table [5] and compared to the 2006 PDG averages © 
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Figure 2: Ratio of preliminary D hadronic branching fractions to the 2004 PDG 
averages (dots). The shaded bars represent the errors in the PDG averages. 
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Figure 3: Afbc distribution for Df — > K + K ir + events. The narrow peak 
at 2.04 GeV/c 2 is due to D s produced directly, while the broad peak between 
2.01 - 2.07 GeV/c 2 is due to D s from D* decay. 



K S K* 



K* K" it* 



K* K" if TC° 
tC jt* it" 
jt* r| 
Jt* t]' 

1 2 3 4 5 6 7 
Branching Fraction (%) 

Figure 4: Preliminary D s branching fractions (dots with error bars) compared 
to the 2006 PDG averages (shaded bars). 

in Fig. |H 

Table 2: Preliminary D s branching fractions with statistical and systematic 
uncertainties. 



Decay 




D+ -> 




1.50 ±0.09 ±0.05 


D+ -> 




5.57 ±0.30 ±0.19 


D+ -» 




5.62 ±0.33 ±0.51 


£>+ -» 


7r + 7r + 7r _ 


1.12 ±0.08 ±0.05 


£>+ -> 




1.47 ±0.12 ±0.14 


£>+ -> 


7T + r/ 


4.02 ±0.27 ±0.30 



The decay Df — > 07r + — » Jf + /f _ 7r + , which is one of the largest and 
easiest to reconstruct, is frequently used as a reference mode to normalize 
other D s decays. However, Dalitz plot analysis of this final state by the E687 
and FOCUS collaborations has revealed significant signal contribution (from 
/°(980) or a°(980)) in the relevant K + K~ mass region. Because of this extra 
signal the 4>ir + branching fraction might be ill-measured depending on the 
specific choice of (mass and helicity angle) cuts. Therefore, we report the 
partial D+ — > K + K~ir + branching fraction (Bam) where the mass of the 
K + K~ system lies within a ±AM (in MeV/c 2 ) mass range around the <f> mass 




(1019.5 MeV/c 2 ), which is more useful from experimental point of view than 
the <pTT + branching fraction. The partial branching fraction with two choices 
of AM are B w = (1.98 ± 0.12 ± 0.09)% and B 20 = (2.25 ± 0.13 ± 0.12)%. 

3 D + and Df leptonic decays and decay constants /d (s) 

In the Standard Model (SM), purely leptonic decays Df^ — > t+vg proceed 
via the annihilation of the constituent quarks into a virtual W + boson. The 
decay width is proportional to the decay constant, /d (s) , which encapsulates 
the strong interaction dynamics in the decay: 

- i + v) = ^m 2 M (l " I W/£ (s) , (2) 

m and M are the mass of the charged lepton and the Dr a \ meson, respectively, 
Gf is the Fermi coupling constant, V c d (V cs ) is the relevant quark mixing 
(CKM) matrix element. 

Knowledge of the decay constants is critical for the extraction of CKM 
matrix elements: e.g. the determination of Vtd and Vt s from measurement of 
BB and B S B S mixing is limited by the uncertainty in the calculation of fs and 
Jb 3 , which currently cannot be measured directly. Experimental measurement 
of the D meson decay constants (/d (s) ) provide an important test of strong in- 
teraction theories and validate the most promising calculations involving lattice 

QCD ©. 

Since the decay width is a function of m 2 (helicity suppression), the decay 
rate to rv is the largest among the three lepton flavors. Although the D + 
(Df) decay rate to [iv is a factor of 2.65 (9.72) smaller in the SM, it is easier 
to measure than the decay to rv because of the presence of extra neutrino(s) 
produced by the subsequent decay of the r. The decay rate to ev is suppressed 
by about five orders of magnitude which is well below the current experimental 
sensitivity. Any deviation from the SM ratios would be an indication of new 

physics 

CLEO has measured the D + —> n + v branching fraction in 281 pb _1 data 

collected at the ^(3770) ®. We have fully reconstructed the D~ decaying to 
sixhadronic final states (K + ir~Tr~ , K + Tr~Tr~ir° , KgW~, Kgir~ir~Tr + , Kgir~Tr°, 
K + K~tt~), which represent more then 35% of all D decays. Candidates are 
selected by requiring AE to be consistent with zero within 2.5<tab, and the 
number of tags in each mode is extracted from a fit to the M^ c distribution. The 
sum of all tags in the range — 2.5cr Mbo < M hc — M D < 2.0oM bc is 158, 354 ±496 
with a background of 30, 677. In the selected events, we search for a single 
additional track consistent with a /i + and calculate the missing mass squared 

MM 2 = (£ boam - E^ + f - (-p D - - p, I+ ) 2 , (3) 
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Figure 5: The MM 2 distribution for D + — ► [i + v candidate events in data. The 
insert shows the region around zero where the arrows indicate the ±2a signal 
region. 



where pjj- is the three-momentum of the fully reconstructed D~. The MM 2 
distribution for the data is shown in Fig. The peak near zero is mostly due 
to D + — > n + v signal, while the peak at 0.25 GeV 2 is from D + — > K°ir + decays 
when a Kl escapes detection. 

The signal region within 2a around zero contains 50 events and the total 
background is estimated to be 2.8 ± 0.31q3 events. After background sub- 
traction and efficiency correction, the measured branching fraction is B(D + — > 



= (4.40 ± 0.66±g:?|) x 10" 4 . The decay constant obtained from Eq. 



H using \V cd \ = 0.2238 ± 0.0029 and the D+ lifetime (1.040 ± 0.007 ps) 
f D = (222.6 ± 16.7±|;1) MeV. 



IS 



We also search for D + 



o + 



v decay by requiring that the extra track is 



consistent with an electron and set a 90% C.L. upper limit of B(D + 
2.4 x 10~ 5 in the absence of any signal. 



< 



and Df 



The branching fraction of Df — ► fx ' v 
measured in 314 pb _1 data collected at e + e~ collision energy near 4170 MeV 
We fully reconstruct one Dj in eight hadronic decay modes (K + K~tt 



is 



TjTT 



K Q S K- , 

7T+7T 7r , K* K*°, rip ). Tags are selected by requiring the 



beam constrained mass to be 2.015 



VP 

< M hc 



< 2.067 GcV/c 2 which is wide 



enough to accepts both direct as well as secondary D s from D* decay. The 
number of tags in each mode is extracted from a fit to the invariant mass 
distribution of the D~ candidates. There is a total of 31, 302±472 reconstructed 
tags within 2.5a (2a for the rjp~ mode) of the D s mass. In contrast to the 
hadronic branching fraction measurement, we select a 7 candidate assumed to 
be the photon from the D* — > jD s decay, and then calculate the recoil mass 
against the Dj tag and the 7: 

MM* 2 = {E cm - E Ds - E 7 ) 2 - (p cm - p Ds - p 7 ) 2 , (4) 

where E cm (p C m) the center of mass energy and momentum of the colliding 
e + e~ beam. Regardless whether the D~ candidate is from the D* decay or 
not the recoil mass should peak at the D s mass. We use kinematic constrains 
to improve the mass resolutions and remove multiple combinations. The recoil 
mass spectrum of each decay mode is fitted individually to extract the number 
of D*D S candidates, which result in a total of 18, 645 ± 426 events within 2.5<r 
interval around the D s mass. The invariant mass and recoil mass distributions 
for Dj — * K + K~tt~ candidates are shown in Fig. O 
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Figure 6: Invariant mass of Dj — > K + K~ir~ tags (left) and the recoil mass 
against the same tag and an additional 7 (right). 

Then we require a single additional track in the event with opposite charge 
to the D s tag and no additional neutral energy cluster with more then 300 
MeV. The missing mass is calculated using the energy and momentum of the 
candidate track (E^, p^): 

MM 2 = (E cm - E Ds -E y - E^) 2 - (p cm - p Ds -p y - p^) 2 . (5) 

We consider three cases depending on whether the additional track is consistent 
with (i) muon (from D s — > /iv), or (ii) pion (from D s — > tv — > TtvD), or (iii) 



electron (from D s — > ev) . The separation between our muon and pion selection 
is not complete: the muon selection is 99% efficient for muons (with a 60% fake 
rate for pions), while the pion selection accepts 40% of pions (with a 1% fake 
rate for muons). The MM 2 distribution for the three cases is shown on Fig. 
[71 The peak around zero in (i) is mostly due to D s — > fiv events. In contrast, 
D s — > tv — > -KVV events has a longer tail on the positive side due to the extra 
neutrino. Therefore, we define three signal regions: (A) —0.05 < MM 2 < 0.05 
GeV 2 in (i) for fiu (92 events); (B) 0.05 < MM 2 < 0.20 GeV 2 in (i) and (C) 
-0.05 < MM 2 < 0.20 GeV 2 in (ii) for ttvv (31 and 25 events, respectively). 
The estimated background from sources other than D s — > \w or irvv decays is 
3.5, 3.5, and 3.7 events, respectively, in the three signal regions. 
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Figure 7: The MM 2 distribution in data when the additional track is consistent 
with muon (i), pion (ii), or electron (iii). 

We calculate three branching fractions: B(D+ — > n + v) = (0.594±0.066± 
0.031)% using signal region (A) only; B cS {D+ -> (i+v) = (0.621 ± 0.058 ± 
0.032)% from combining all three signal regions (A) + (B)+(C); and B(D+ — > 
t+v) = (8.0 ± 1.3 ± 0.4)% from the combined tv regions (B)+(C). In the 
first two cases, the D s — » tv contribution is subtracted assuming the relative 



decay rate between \xv and rv final states is equal to the SM expectation and 
using B(t -» ttv) = (10.90 ± 0.07)%. We also set a 90% C.L. upper limit on 
B(D+ -> a+v) < 1.3 x 10~ 4 . 

The decay constant is calculated from the most precise branching fraction 
(B cS ) using Eq. EH with \V CS \ = 0.9730 and the D s life time of (500 ± 7) x 10~ 15 
s: f Da = (270 ± 13 ± 7) MeV. 

We also measure Dj — > t + v with a different technique utilizing the r + — > 
e + vv decay with a total product branching fraction of about 1.3%. In this case, 
we fully reconstruct the Dj candidate in the event and require an additional 
track consistent with an e + but do not attempt to find the 7 from the D* decay. 
Events with additional tracks and more than 400 MeV total neutral energy in 
the calorimeter are vetoed (the typical energy of the 7 or ir° from D* decay is 
around 150 MeV). After analyzing 195 pb~* subsample of our data we obtain 
a preliminary branching fraction B(D+ — » t + v) = (6.29 ± 0.78 ± 0.52)% and 
decay constant f Dg = (278 ± 17 ± 12) MeV. 

The weighted average of these two results is fr> 3 = (273 ± 10 ± 5) MeV. 
Combined with our published fo value we find a ratio JdJ Id = 1.22 ± 0.09 ± 
0.03. The measured decay constants are consistent with most theoretical mod- 
els. In particular, recent unquenched Lattice QCD calculations ® yield fu = 
(201±3±17) MeV, f Ds = (249±3±16) MeV, and f D Jf D = 1.24±0.01 ±0.07. 

4 Measurement of D° mass 

Precise knowledge of the D° mass is not only important for its own sake but 
it can also help with the interpretation of the V(3872) state. Because of the 
proximity of the X mass (3871.2 ± 0.5 MeV/c 2 ) to M(D°) + M(D*°), one 
theoretical suggestion is that the X(3872) is a bound state of D° and D*° 

mesons However, it is necessary to measure the D° mass with better 

precision than the current PDG average of 1864.1 ± 1.0 MeV/c 2 ^ in order to 
reach a firm conclusion. 

CLEO has measured the D° mass ELll hi e + e - ^ ^(3770) D°D° 
events using the decay D° — * Kg(j) followed by Kg — >_7r + 7r _ and (f> — > K + K~ . 
In order to obtain a clean sample of signal events, the D° has been reconstructed 
using the same tagging technique described in section^ imposing loose require- 
ments on AE and M^ c of the candidates. The D° — > Kg4> decay was selected 
because the final state pions and kaons have small momenta, and therefore the 
uncertainty in their measurements makes small contribution to the final result. 
In addition, the mass of the Kg candidates can be kinematically constrained 
to its well known value. 

Pions from the K s are required to originate from a displaced vertex and 
have a M(tt + it~) invariant mass in the range 497.7 ± 12.0 MeV/c 2 before 
the mass-constrained kinematic fit. The (j) candidates are accepted with a 



M{K + K~) invariant mass of 1019.46 ± 15 MeV/c 2 . The mass distribution of 
the D° candidates in 281 pb _1 data is shown in Fig. [H] A likelihood fit using 
a Gaussian peak and a constant background yields 319 ± 18 D° events and 
a D° mass of 1864.847 ± 0.150 MeV/c 2 with a mass resolution of 2.52 ± 0.12 
MeV/c 2 (the errors are statistical only). The total systematic error on the 
mass measurement (0.095 MeV/c 2 ) is dominated by uncertainty in detector 
calibration, which is studied using the K$ mass in inclusive D — > KgX decays 
and the ip(2S) mass in exclusive ip(2S) — > tt + tt~ J/tp(J/tp — > /j + aO events. 
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Figure 8: The invariant mass of D° — ► K®(j> candidates in data. 

Our final D° mass with statistical and systematic uncertainties is 

M(D°) = 1864.847 ± 0.150 ± 0.095 MeV/c 2 . (6) 

This gives M(D°D*°) = 2M(D°) + AM D >o_ D o = 3871.81 ± 0.36 MeV/c 2 , 
and leads to a binding energy of the X(3872) as a D°D*° molecule: AEb = 
M(D°D*°) - M(X) = +0.6 ± 0.6 MeV/c 2 . The error in the binding energy is 
now dominated by the uncertainty in the mass of the JT(3872). 

5 Study of Xcj —> h + h~h° decays 

In contrast to the 1 members of the charmonium states (J/tp, ip(2S)), the 
decays of the Xcj (J — 0,1,2) states are not well studied. The different decay 
mechanism of these states (dominated by annihilation into two (virtual) gluons 
and contribution from the color-octet mechanism) might provide complimen- 
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tary information on light hadron spectroscopy and possible glueball dynamics 

At CLEO, the Xcj states are produced in radiative decays of the ip(2S) 
and we study their decays to eight selected three-body hadronic modes: TT + TT~r], 
K + K~r), ppn, ir + Tr~T]' , K + K~ir°, ppir , ir + K~Ks, and K + pA. We have mea- 
sured branching fractions or set upper limits for the first time in most cases 

using about 3 million ip(2S) decays . As an example, Fig. [H] illustrates the 
invariant mass distribution for two of the hadronic final states. 
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Figure 9: The invariant mass distribution for Xc.j — * K + pA (left) and Xcj 
Tr + Tr~rj (right) candidate events in data. 



We perform a Dalitz-plot analysis of the decays with the highest statis- 
tics, Xd -> 7T+7r-?7 (228 events), K + K~n° (137 events), and tt + K~ K% (234 
events) in order to study the two-body resonant substructure. We use a sim- 
plified model with non-interfering resonances, which is adequate to show the 
largest contributions in our small sample. Fig. [TO] shows the Dalitz plot and 
three projections for Xci — * 7r + 7r~?7 and the result of the fit. There are clear 
contributions from ao(980) ± 7r =F and /2(1270)?7 intermediate states, and a sig- 
nificant accumulation of events at low tt + tt~ mass which can be described by 
an S-wave (tr) resonance. This mode might offer the best measurement of the 
ao(980) parameters with higher statistics. The decays Xci — * K + K~ir° and 
tt + K ~Kg are analyzed simultaneously taking advantage of isospin symmetry. 
We observe contributions from K*(S92)K, ^(1430)^ , a (980)7r intermediate 
states. It is not clear whether the K*(1A30) is Kq or K%, and other Ktt and 
KK resonances can contribute. Addition of kK or non-resonant component 
does not improve the fit and the significance of their contribution remains under 
3 standard deviation. 

More data is required to do a complete partial-wave analysis taking into 




Figure 10: Dalitz plot and projections of Xci ^ vr + 7r 77 decay. 



account the Xci polarization properly and including interference among the 
resonances. 



6 Summary 

I have reported mostly preliminary results for hadronic and purely leptonic 
decays of D and D s mesons from the CLEO-c experiment. These results rep- 
resent substantial improvement over previous measurements. The D + and D° 
hadronic branching fractions are limited by systematic uncertainties of up to 
3%. The D s hadronic branching fractions are measured with relative uncer- 
tainties between 6 — 12%, which are dominated by statistics. The measurement 



of the decay constants fo (/r> (s) ) from purely leptonic decays are also statistics 
limited with a total relative uncertainty of 8% (4%). 

I have also presented the most precise measurement of the D° mass, and 
a study of three-body hadronic decays of the \ c j states. 

Precision of these and other measurements will improve in the near future 
with more data on the way. CLEO-c has already collected an additional 8 times 
more data on the tp(2S), and we plan to increase the DD and D^Df data 
samples by a factor of 2-3 before data taking ends in April 2008. 
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